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Exact analytic expression is derived for the matrix elements of the Coulomb interaction in two 
dimensions in the form of a closed finite sum expression. The orthonormal complete set of eigen- 
functions of the harmonic oscillator | n x n y ) is used as the basis for spanning real space. Several 
recurrence relations have been found in order to simplify the task of calculating the usually vast 
amount of elements required for any computer simulation. 
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The choice of the anisotropic orthonormal complete set of eigenfunctions of the harmonic oscillator is due to the 
ubiquity and special symmetric properties of this particular problem in the field of quantum mechanics. The projection 
of this basis into real space is 
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Without any loss of generality we can assume that 
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O ■ ax> a v ( 3 ) 



Just the case a x = a y will fall apart from our supposition. There is no problem in treating this case apart, as it is 
even simpler to calculate and can be derived from our results. The eccentricity parameter 7 is defined as 
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J* < 7 < I- (5) 

> ■ 

It is of great interest to express the wavefunctions and Coulomb potential as their two-dimensional Fourier trans- 
forms 

Mr) = ^ f MQ)e- l "- r dq, (6) 
k> , 2ir J 

V(r) = -L J V^e-^dq. (7) 
Under this notation, the matrix elements take a different expression 

= ^ / - + Q)V(q) d qi dq 2 dq, (8) 

where Xi represents the set of quantum numbers {n x , n y } for the i-th particle. 

<t = ^ / {q)D% (q)V(q) dq, (9) 

where 

C X x ,{q) = J cj>l{k)<t> x ,{k-q)dk= [rx(r)M(r)e- iq - r dr, (10) 

Dy(q) = [ <t>l(k)<t>X'(k+q)dk= J ^xi^y (r)e^ r dr 

= Cl,{-q). (11) 



The calculation of C\\>(q) is quite straightforward. It can be separated in both two dimensions, as they are (almost) 
symmetrical. The calculation proceeds as follows (normalization constants omitted for clarity) 

/oo />oo 
e- x2 / a2 *H nl Jx/a x )H nt (x/a x )e-^ x dx = a x / e~ x ' 2 H n i (x')H nt (x')e~^ a * x ' dx' 
-OO J —GO 

/OO 
e- x ' 2 H ni (x')H ni (x')e-^ x ' dx', (12) 
-OO 

where x' — a x x and q' x — a x q x . This equation can be seen as the Fourier transform of the product of two functions: 

f [/</<]> ( 13 ) 

being 

fn(x) = e- x2 / 2 H n (x). (14) 

Making use of the properties of the Fourier transform, in particular the one that states that the transform of the 
product of two functions is the convolution of their transforms, i.e. 

In our particular case, the Fourier transform of the function f n (x) has a very simple form as it is an eigenfunction of 
this operator with associated eigenvalue i n : 

T[f n {x)](q x ) = l n f n {q x ). (16) 
It only remains now to calculate the convolution product, defined as 

/OO yOO 
f(r)g(t-T)dr= / g(r)f(t - r) dr. (17) 
-00 J — OO 

In our case, we have to calculate the following integral 

/OO 
e- x ' 2 ' 2 H nl {x')e-^- x ^' 2 H nt { q ' x - x')dx'. (18) 
-oo 

Let us first do the following change of variables: 



x 1 = u + q' x /2 q' x = 2x' - 2u 

dx' = du (— oo, oo) — > (— oo, oo) 



(19) 



Thus, choosing n x \_ = m&x(n x ,n x ) and n x A _ — m.m(n xl n x ) , the last integral is transformed to 

/oo 
e -(u +q 'j2 f/2 Hnlt {u + q 'J 2 )e-^-^ 2 l 2 H nl% (q' x u q' x /2) du 
-oo 

/OO 
e- u2 H <t (u + q' x /2)(-lT^H nlx (u - q' x /2) du, (20) 
-oo 

which according to Gradshteyn & Ryzhik Eq. 7.377 in page 797 

/oo 
e- x2 H m (x + y)H n (x + z) dx = 2 n ^' 2 m\z n ~ m L n ^ m {~2yz) [m<n], (21) 
-OO 

the result of integral (20) is 
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9 n" 1/2 14 , ( axq x \ n * + - n *- ,n&-»i- 2 2 



3 



Collecting all terms together, we can write the final result for C(q): 

1/2 



C <nl^<ly) ( axay7r ) ( 2 »in£! 2 n ln\\ 1<n%\ 2<n y \ 



a x a y 2 n *+ 2"i+ \n y t ! (- 1 )""+ +""+ 

e-^^r-^r" + " B "^- ,,lt -(^/2) 



I$r ni - {alql/2)L<i-<- (a 2 y q 2 y /2). (23) 



In the same fashion we calculate the function D: 



_ 2 n- n 23_, 2 n- „23_,\ V2 



Dft(q x ,q v ) = (-1)"*++S+ " v 



23 23 23 23 

^3 + -"M«^/2)^+-"M^/2)- (24) 

as- y- 

The Fourier transform of the Coulombic potential is 

V(q) = / V(r)J (rq)rdr = f- = £— / e -(«2+<£)« a du (25) 
Jo 9 "v/Wo 

2/' 00 /22 l 222/2\2/2 

= £ / e -K9.K«„« 9 /« 9 )" /« I( i M (26) 

= £a x 4= / e-^ 2+ «i 2 ^ 2 )" 2 dn (27) 
V* Jo 

V(r) = i (28) 
Where we have substituted the dimensionless momentum variables, 

q' x = a x qx, (29) 
q' y = a y q y . (30) 

The integration can be now performed. Making use of the expansion of the Laguerre polynomials as 

k=0 v 7 

let us focus on the q x and q y dependent part. Let us also forget about any multiplicative constants and summation 
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indices, we can write the matrix element as: 

(_l)s*+s»+™ii+™"+« 2 J 3 _+™ 2 , ; L 
f 2 



"x- \'"x 



1 \ s x -\-s y J rk x -\-k x -\-ky J rk /■ oo />og 

2 J X loo ^ 

oc J „/ 

r^. e -(« 2 /7 2 + l/2) 9 ; 2 g ' (2s y +2fe y +2^) (- 32 ) 
-oo a V 

The last three integrals can be solved quite easily. First, the ones depending in q' x and lead to: 

/OO -I 
^ ^+1/2)^, (2 M+ 2 K) ^ = r+1/ 2 T {a + 1/2) _ _ _]__ (33) 

/OO -i 
^ ^v^wy^.*^**;) d(?; = 2&+ i/ 2r (6 + 1/2) (i + 2m2 ^ 2)&+1/2 (34) 

where a = s x + fcj, + k' x and b = s y + k y + k' y . We have already included the condition for the integral to be non-zero: 

2s x = \nl-ni\ + \nl-nl\=n x %-n x t+nf + -nf_, (35) 
2s y = \nl-n*\ + \nl-n 3 y \ = n\\ - n\t + nf + - nf_, (36) 
s x ,s y e N. (37) 

Now, the u-integration can be performed. 

1 1 .1/2 \~ 1/2 /l . 1 



o (l + 2u 2 / 7 2 ) b + 1 / 2 (i + 2?i2 )a+i/ 2 du 2U 2 / B U' a + 6+ 2 

1 1 



x 2 f\ ( a+ l l;a + 6+l,l- 7 2 ) (38) 



which was solved with the help of Gradshteyn & Ryzhik Integral 3.259 3, page 326 (check also the errata) 

/ x A - x (l + ax p )-»(l + (3x p )~ v dx = -a~ x / p B (-^ + v --\ 2 F X ( v, -; p + v; 1 - ^] 
Jo P \P P) V P a J 



(39) 



[|arga|<7r, |arg/3|<7r, p > 0, < Re A < 2Rc(^ + v)] 
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Grouping all terms together, we end up with: 

y n * n l n2 x n y _ J_/_iyi£'_-r». " " . 



1 , ^23_ 23 /2< 4 +nii!2"" nii!2"" n 23 ^ 2"?+ n 23 ! \ 1/2 



"x '"v '"a: '"y 



2tt v ^ \^ 4 + ! 2< 4 - n&! 2»" ^ 3 + ! 2<- nf + \ T f- 

(-iy 
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- 14 n 23 

(-l) fe * / < + \ ^ (-1)*, / n 14 



V v ~ ' y+ V ZZ 

Z^ kJ In 14 - fc. / Zw y 



k Vj,- '"!// fe , ■»!/■ Vj,- '"3/ 



X " " — — 2^+ s y+ k -+ k y+ k ^+ k 'y +1 r( Sx + k x + k , x + l/2)T(s v + k y + k' y + l/2) 
y z y &y 

K"2") B (i Sx + ^ + ^ + ^ + fc - + ^ + 

2F1 + k x + k' x + 1 1; s x + s a + fc, + fc y + fc; + fc; + 1, 1 - 7 2 ) • 



(40) 



After some compactation we get the final result: 



1 f 232 3 /n 14 I n 14 1 n 23 I n 23 '\ 1/2 



a, V^ 3 ' \nl\\nl\\nf + \nf + \) 

rc 14 n 23 

B ^— , s x + Sj, + k x + fcj, + k x + fc^ + — ^ 

2*1 (s, + + + ^, % + s„ + fc, + ky + k' x + k' y + 1, 1 - 7 2 ^ . 



(41) 



RECURRENCE 



In order to optimize the computing times, let us try to find some recurrence relations for the kernel of the summation. 

Vij = B(±,i + j + ±) 2 F 1 (i+ §, i;z + j + l;l- 7 2 ) 

= B i+j F iti+j . (42) 



The Beta function behaves as follows 



B = tt, (43) 
B k = ^^B fc _i. (44) 
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The Gauss hypergeometric function does it in a more complex way 

a(z - l) 2 Fi(a + 1, b; c; z) = (a - c)2-Fi(a - 1, b; c; z) + (c - 2a + (a - 6)^)2-^1(0, b; c; z) (45) 
= (c - l)c(«- 1)2^1(0, 6, c-l,«) 

+c(c - 1 + (a + b - 2c + l)z) 2 Fi (a, 6, c, z) 

+(a- c)(b- c)z 2 F 1 (a,b,c+ l,z) (46) 
, , \ (c-l)(2-c-(a + b-2c + 3)z) „ , , 
»*(°.^*) ^ ( fl _ c + 1 ) (6 _ c+1)z — 2-Pi(a.&;c- 1;^) 

+ / C - 1 | ( ;- 2)(1 -^ ) 2 F 1 (a,6; C -2 ; z), (47) 
(a — c + l)(o - c + l)z 

„, . c- 2a + 2+ (a - - l)z 

2 Fi(a,6;c;z) = — 2 Fi(a - 1, b; c; 2) 

(a - l)(z - 1) 

2 Fi(a-2,6;c;z). (48) 



(a-l)(*-l) 

In our case a — k+\, b=^ and c = / + 1. Therefore if F fe; = 2 Fi(fc + |, |; / + 1; z): 

_ l{l-l-(2 + k-2l)z) T? 

~ a-i)a + k-i)z 



-^- 2 , (49) 



(i-Z)(i + fc-Z> 
2-2fc + / + (fc-l)z F 
(fc- 1) 

where z = 1 — -f 2 . Now combining expressions (44), (49) and (50), it is possible to find a recursive relation for the 
rows and columns of Vij. For the columns: 

Vij = Bj+ji^j+j 

( (* + j)(l-*-J-(2-i-2j)(l-7 2 )) „ 

(i+j-l)(i+j)7 2 F 
_ l-i-]-(2-i-2 ] )(l- 1 2 ) (»+J-l)7 2 B p 

D l+ j-i-tiA +: j-i + — — - — Uj + j_i^i ii+ j_ 2 

(z+. 7 -i) 7 2 » + j - § 

w *,7"-i + ttt: ~ ; — : — rt>i 





-I 2 ) 


1 - i ~ 3 ~ (2 - i 


-2j)(l-7 2 ) 




_ 7 2) 


1 - * - i - (2 - i 


-2j)(l-7 2 ) 



(i-i)(i-7 2 )* + i-i 



(i + j — -)t 2 

+ 1U1 2^ -2' ( 51 ) 



U-i)(l-7 2 ) ^ (i-|)(l-7 2 ) 



And for the rows: 



/ l+j + (l-i) 7 2 (j+j) 

- Hi+j I ( 1 _ i) 7 2 + (i _ 1 ) 7 2 **-2.*+j 



l+,7 + (l-z ) 7 2 , 3 + - 



Vi-l,j+l + T T r-^V i -2j+2- (52) 



,2 



i) 7 2 ' (i - |) 7 



The diagram for the column iteration (with fixed row). 



, - 2 ■ Vj J 1 — Vjj 



And for row iteration. 



v j 



Vi-2,j+2 




< > %' • • 

The seed for the iteration is: 



Vij 





1 2 





2X(z) 




1 


\{K{z)-E{z)) 


4, (2(z - l)K(z) - (z - 2)E{z)) 


2 







NOTES 

The Gamma function for half-integer values is defined as: 

_ / n (2n-l)!! r- 

r (^+l)= 2n 
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